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I. INTRODUCTION 

High energy exclusive leptoproduction processes have 
been drawing increasing attention after a long hiatus 
since they were first proposed as direct probes of par- 
tonic structure. The first exclusive electron proton scat- 
tering experiments were conducted at both DESY (HI, 
ZEUS and HERMES) and Jefferson Lab. A new forth- 
coming dedicated set of experiments are currently be- 
ing performed and planned at both Jefferson Lab 1 and 
CERN (Compass) [2]. The possibility of using neutrino 
beams to study hard exclusive reactions also concretely 
exists within e.g. the Minerva experiment at Fcrmilab 

The interest in deeply virtual exclusive processes orig- 
inates from the realization first discussed in Refs.[l]-[5] 
that QCD factorization theorems similar to the inclusive 
DIS case can be proven. CoUinear factorization theorems 
have in fact so far been established for Deeply Virtual 
Compton Scattering (DVCS), involving initially trans- 
verse photons, and for Deeply Virtual Meson Production 
(DVMP), with initial longitudinal photons New de- 
velopments in QCD factorization are also rapidly evolv- 
ing [S]. The leading order diagrams describing the am- 
plitude for the scattering process are shown in Figure [T] 
A phenomenology ensues similar to the one extensively 
developed for inclusive scattering, an important differ- 
ence being that exclusive reactions provide us with addi- 
tional kinematical dependence on the momentum trans- 
fer squared between the initial and final proton, t, and on 
its Light Cone (LC) component, C,. The new t-channel 
variables allow us to pin down in principle the depen- 
dence of the parton distributions on spatial d.o.f. through 
Fourier transforms of the Generalized Parton Distribu- 
tions (GPDs) [9 . The latter enter the description of the 
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FIG. 1: Leading order amplitude for the DVCS process, 7* -I- 
P 7 + P' 



soft matrix elements for DVCS, DVMP and related pro- 
cesses. The kinematical variables external to the partonic 
loop in Fig{l^ are the initial photon's virtuality, , the 
skewness, C = {Aq)/{Pq) ~ Q^/2{Pq) = xbj, t = A^, 
with A = P - P', P (P') being the initial (final) photon 
momentum. (C,i, Q^) define a set of independent invari- 
ants. In the factorized approach one defines also internal 
loop variables, X = {kq)/[Pq) - the parton's momentum 
fraction - and fc^ - the intrinsic transverse momentum. 

At high momentum transfer, the amplitude for DVCS 
can be written schematically as 

q 

(1) 

where we considered for ease of presentation only the 
GPD H (detailed expressions will be given in what fol- 
lows). The analog of the Compton Form Factor (CFF) 
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with one virtual photon is 



+1 



-i+C 



X-C + ie X-ie 



(2) 



One has therefore that both the imaginary and real parts 
of the amplitude, namely 



P.V. j ^dXH,{XX,t) 



X-C + ie X 



(3a) 



(3b) 



enter the description the DVCS reaction. Information on 
the partonic distributions which is contained Hq{X, C,t), 
needs to be extracted from these observables. This is an 
important difference with DIS where, because of the opti- 
cal theorem, the cross section by definition measures the 
imaginary part of the forward amplitude. The DIS cross 
section is therefore directly proportional to linear com- 
binations of the soft matrix elements, or Parton Distri- 
butions Functions (PDFs) convoluted with appropriate 
Wilson coefficient functions. On the contrary, in DVCS, 
DVMP, and related processes one needs to disentangle 
both the real and imaginary contributions of the CFFs 
defined in Eqs.( |3a|3b| [10^. 

It was recently suggested to use dispersion relations 
in order to relate the real and imaginary parts of CFFs. 
However, as we noted in [TT], dispersion relations do not 
apply straightforwardly because of the appearance of t- 
dependent physical thresholds. As a result, we reiterate 
that both the real and imaginary parts need to be ex- 
tracted separately from experiment, at variance with the 
simplification suggested e.g. in Refs.|121 fT3] . 

On one hand, the type of information we wish to obtain 
from high energy exclusive experiments is a sufhciently 
large range of GPD values in {(,t,Q^) that would en- 
able us to reconstruct the partonic spatial distributions 
of the nucleon from a Fourier transformation in A . This 
would allow us both to explore the holographic principle 
for the nucleon, and to connect to Transverse Momentum 
Distributions (TMDs). On the other hand, it is impor- 
tant to have access to the spin flip GPD, E, which is es- 
sential for determining the Orbital Angular Momentum 
(0AM) contribution to the spin sum rule . 

The question of whether the various GPDs can be 
extracted reliably from current experiments has been 
raised, given the complications inherent both in their 
convolution form, and in their complex multi-variable 
analysis (see e.g. [Ml [15]). A pragmatic response was 
given in [161 117] where an assessment was made of which 
GPDs can be extracted using the present body of data 
from Jefferson Lab and Hermes. In particular, it was 



concluded that the only CFFs that are presently con- 
strained by experiments are HeH and IrnH, with rather 
large errors, up to 30%. A global fit using the dual model 
of Ref. [13], valid mostly at low Bjorken x was also con- 
ducted in |15j . However, these approaches raise many 
concerns. In particular, can the "dual model" used in 
the fits accommodate all of the data with the given num- 
ber of parameters? The problem is critical, in particular, 
for both higher C values and for the real CFFs. Further- 
more, the analysis of [17] does obtain model independent 
extractions of CFFs at the expense of not allowing for 
extrapolations to kinematical domains beyond the very 
sparse data sets. 

The goal of our fit is to extract the GPDs from a variety 
of experiments under the following basic assumptions: 

i) QCD Factorization is working, namely the soft and 
hard parts are separated as shown in Figjl] 

ii) the GPDs contributing to DVCS are evaluated at the 
lowest order in QCD. 

This situation is somewhat similar to the extraction 
of PDFs from structure functions at NLO, where the 
PDFs are convoluted with the NLO Wilson coefficient 
functions. The convolution is neither a substantial or 
conceptual obstacle so long as one is providing an appro- 
priate initial functional shape. The strategy we propose 
here provides a parametric form of the chiral even GPDs, 
H, E, H, E that is valid in the multi-GeV, intermediate 
XBj region accessible at Jefferson Lab and COMPASS. 

We suggest the idea that for extracting GPDs from ex- 
periment a progressive/recursive fit should be used rather 
than a global fit. In our fitting procedure constraints 
are applied sequentially, the final result being updated 
upon including each new constraint. In a nutshell, in a 
first step we provide a flexible form that includes all con- 
straints from inclusive data - DIS structure functions and 
elastic electroweak form factors. We subsequently eval- 
uate the impact of presently available DVCS data from 
both Jefferson Lab and HERMES. The data set used in 
our analysis is consistent with the one from Ref. (16j . The 
parametric form is based on a "Regge improved" diquark 
model that because of its similarities and possible connec- 
tions with the dual model |18j , we call the hybrid model. 

Our approach however makes two important distinc- 
tions: i) we attack the GPD parametrization issue 
from the bottom-up perspective. We adopt a flexible 
parametrization that is consistent with theoretical con- 
straints imposed numerically, and let the experimental 
data guide the shape of the parametrization as closely 
as possible. In this procedure, experimental evidence 
is used to constrain the various theoretical aspects of 
the GPDs behavior, eventually giving rise to a complete 
model; ii) our model differs from some of the lore on the 
partonic interpretation of GPDs in the ERBL region. In 
Ref. [in] we in fact pointed out that the ERBL region, or 
the region with X < (, cannot be described in terms of 
a quark anti-quark pair emerging from the nucleon be- 
cause of the presence of semi-disconnected, unphysical, 
diagrams associated to this configuration (Fig{T]D) . While 



3 



casting a doubt on any partonic picture in the ERBL re- 
gion, we also suggested that multi-parton distributions 
may restore the connectedness, and consequently the par- 
tonic interpretation of the DVCS graphs. In this paper 
we therefore adopted, as a practical scheme, the hybrid 
model in the DGLAP {X > () region, and a minimal 
model that is consistent with the properties of continuity 
at X = C, polynomiality, and crossing symmetry for the 
ERBL region. 

Details on the model are given in Secti on [iTl The new 
fitting procedure is described in Section |III| In Section 
|IV| we discuss the implementation of world DVCS data. 
Finally, in Section[V]we draw our conclusions and outline 
future work. 



II. COVARIANT FORMULATION AND 
SYMMETRIES 



respect to /i o z/, that will enter the unpolarized and 




We begin by describing the connection between the 
Dirac basis formulation of the correlation function and 
the helicity amplitudes formalism in DVCS. Some of this 
formalism was outlined in Refs.pUl [2T] . We, however, 
present here the formal details that will be important for 
the extraction of observables in Section ITVI 



A. Formalism 

The factorization theorem for hard scattering processes 
allows us to separate the hard scattering between the el- 
ementary constituents, which is calculated using pertur- 
bative QCD, from the soft hadronic matrix element, Ai, 
as 



+ 



rTr 



(27r)4 " [V (fc + g)2 +ie 



Ai{k, P, A) is the off-forward correlation function: 



(4) 



Mtf{k,P,A) 



d%e^'y{p',A'\^J0)My)\P,^)A5) 



where we have written out explicitly both the Dirac in- 
dices i, j, and the target's spins A, A'. By projecting out 
the dominant contribution in the Bjorken limit (Q^ — cx), 
XBj = /IMv C and t fixed), which corresponds to 
transverse virtual photon polarization, one obtains 



FIG. 2: Left: Amplitude for DVCS at leading order in 
The light cone components of the momenta for the active 
quarks and nucleons are explicitly written; Right: Time 
ordered diagrams for DVCS: (a) dominant contribution in 
V > ^ region; (b) a qq pair is first produced from the nu- 
cleon and subsequently interacts with the photons. This pro- 
cess dominates the X < C, region. The crossed-terms where 
two of the particles in the same class are switched, are not 
shown in the figure. 



longitudinally polarized scattering, respectively, as will 
be clarified in what follows. 

A possible kinematical choice is the one where the 
struck quark's light cone (LC) longitudinal component 
is fc"*" = VP"*", the momentum transfer, is decom- 
posed into a LC longitudinal component, A+ = CP+ 
and a transverse component, Aj_, such that the invari- 

A^ 

ant t = A^ reads: t = -C'^M'^/{1 - Q) - ^ ^^^^^ 

these variables one can perform an integration over the 
quark loop momentum namely d^k = dk~^dk~d^k± = 
P^dXdk~d^k±, obtaining the following expressions for 
the CFFs, 



dX 



1 



1 



-1+c v^-C- 

X Fi,,{X,C,t), 



X-ie 



(7) 
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where gi^ 



an Su p ^ 



(6) 



= 9 — P n ~ P « , "^T 

e p and n being unit light cone 
vectors. The labels S and A refer to the symmetric and 
antisymmetric components of the hadronic tensor with 



where 



dX - 



i+C 



X 



1 1 

Q + ie ^ 



(8) 



4 



2P 



2M 



(9) 



-Pa,a'(-''^' C, t) - 

2P 



U{P', A') ( ^+j,H{X, C, t) + 15^E{X, C, t) ) U{P, A) 



(10) 



The GPDs, H, E, H, E introduced in the equations in- 
clude the integration over dk-(Pkr, = (P++P+')/2, 
and we did not write exphcitly, for simpHcity, both the 
label for the different quark components, and the depen- 
dence on the scale ol the process, Q^. 

Eqs.(|9]) and ( 10 ) define the basic Dirac structure for the 
chiral-even sector at leading order in 1/Q. The connec- 
tion with the helicity formalism is made by introducing 
the helicity amplitudes for DVCS, 



A^.A:A' .A' 



*A' 



AA''=i^ 



(11) 



where e^, are the photon polarization vectors, (A-y,A) 
refer to the initial (virtual) photon and proton helicities, 
and (A^, A') to the final ones. The following decomposi- 
tion of /a^,A;A' ,A' ED can be made 



/a„A;A;.A' - ^ 5a,V "'(^'C,i;0')<^^A',A';A,A(X,C,t), 

(12) 



where g^^ \, ^ describes the partonic subprocess 7* + g — > 
7 -|- g, i.e. the scattering of a transverse virtual photon 
from a quark with polarization A; Aa',a';A,a is the quark- 
proton helicity amplitude, and the convolution integral 

is given by (X) — > ^^^^i dX. In the Bjorken limit g^^'y^'^ 
reads 



a^.a' 



A^ * A' 



* A' A 



(13) 



with s = (fc + g)2 w g2(X - C)/C and u = (fc' - qf « 
Q^X/C, and q- « {Pq)/P+ = QViKP^)- 

For DVCS one can consider either the sum over the 
transverse helicities of the initial photon, i.e. we take 
it as unpolarized, or the difference of the helicities in 
which case the transverse photon is polarized. The out- 
going photon is on-shell, thus purely transverse or helicity 
±1. So the leading incoming virtual photon will have the 
same helicity in the collinear limit. Only g^'''\.{= gZ'Z 
via Parity conservation) for the direct s pole term or 
g+'+i— 9-]-) for the crossed u pole term will survive. 
For either allowed combination, q^j 



9tX 



9. 



9i+=9tt-9++, 



9+ 



+ ±g^^ = ^X{X - C) f — 4 



± 



X -C + ie X-ie 



(14) 



where we use the S/ A labels for the sum/difference be- 
tween the positive and negative polarized photons for 
the overall process representing the sum/difference be- 
tween the quark states' helicities that arises as the quarks 
emerge from the nucleons. 

The quark helicity or chirality is conserved in this hard 
subprocess for DVCS. Hence the Aa',A';A,a will be chiral 



even. Eq.(14| is the Wilson coefficient from Ref. ^ times 



a kinematical factor that will cancel out when multiplied 
by the soft part as described below. 



The convolution in Eq.( 12 ) yields the following decom- 



position of the transverse photon helicity amplitudes 



J++ 

fl- 
ft- 



/++.++ + /-+,-+ 

5++® (^++,++ + ^-+,-+) 

/++,++ ^ ,/ — h,-+ 

9++ ® - A_+,_+) 

/++.+- + ./-+.- 

.9++® (^-+.++ + ^++,-+) 

,/++.+- ^ /-+. — 

g++ ® (^-+,++ - ^++,-+) 



(15a) 
(15b) 
(15c) 
(15d) 



where, because of parity conservation, A ^ = 

^-+,-+ = A — ,+- = -^*++.-+' and 



++' 



and By calculating explicitly the matrix elements in Eqs. ([t]) 
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and ([S]), using the relations below 



^ -U{P\M)^+U{P,k) - ^^^<5a,a' 



2P 



^^l/(P',A')- 
2P 



l-C/2 
-A^Z/(P,A) = 



(1 - C/2) 
1 



-AAi - 
2M 



5k, 



(16a) 



(16b) 



:C/(P', A')7+75t/(i^, A) = A^^^5a.a<16c) 
2P' ^-Cl^ 



'^^UiP',A')^,^U{P,A) 



c 



2(1 -c/2) 



one obtains the various helicity amplitudes written in 
terms of the following combinations of CFFs for the sym- 
metric part, 



fl- 



l-C/2"- ' (1-C/2)VW 
1 1 Ai+iA2^ 



^/r^ 1 - C/2 2M 



and 



-CV4 



l-C/2 (l-C/2)yW 
C 1 Ai+iA2 



1 - c/2 2M 



f (17a) 
(17b) 



£ (18a) 



(18b) 



for the anti-symmetric component. 

A similar formalism was presented in |20) where, how- 
ever, the helicity amplitudes were identified with the 
combinations, 



/++ + /++ ! f+4 



rA fS I fA fS _ rA 

j++tJ-\ — tJA — 'in — J A — ' 



Here, differently from [20], we distinguish between the 
two possible circular polarizations for the transverse vir- 
tual photon, that generate the S and A components. 
Such components are written out explicitly throughout 
this paper. 



B. The Hybrid Model 



We evaluate the quark-parton helicity amplitudes in 



Eq.( 12 1 using a covariant model. The simplest realization 



of the covariant formalism is the quark-diquark model in 
which the initial proton dissociates into a quark and a 
recoiling fixed mass system with quantum numbers of a 
diquark (Fig.l). The covariant model can be made more 
general by letting the mass of the diquark system vary ac- 
cording to a spectral distribution. Extending the diquark 
mass values generalization corresponds to "reggeizing" 



the covariant model since the spectral distribution can 
then reproduce the Regge behavior which is necessary 
to describe the low X behavior (a likewise scenario was 
considered in the pioneering work of Ref . [22^ ) . Keeping 
this in mind, in this paper we will include a Regge term 
multiplicatively, as shown later on. Because we intro- 
duce similarities, or open possible connections with the 
dual model of Ref. 18 , we denote this model of GPDs 
the hybrid model. 

For reasons explained in the Introduction, we adopt the 
diquark model only for the DGLAP region, where X > 
The DGLAP region can be considered a direct extension 
of the parton model, where the struck quark with ini- 
tial longitudinal momentum fraction X, is reinserted in 
the proton target after reducing its momentum fraction 
to X — C, C being the fraction transferred in the exclu- 
sive reaction. In the DGLAP region the initial and final 
quarks are both off-shell, while the diquark intermediate 
state is on mass shell. The soft part, Aa',a';A,a is de- 
scribed described in terms of GPDs. A is given by an in- 
tegral over the k~ and k± variables (see appendix for de- 
tailed expressions). In the DGLAP region the three soft 
propagators corresponding to the quark that is emitted 
(/c^ — rn^)~^, the quark that is reabsorbed (fc'^ — to^)~^, 
and to the intermediate diquark system {P^ — A/|-)~^, 
have poles that lie respectively on the negative imaginary 
axis (fc, k'), and on the positive axis (Px)- Therefore 
one closes the integration contour on the positive side, 
and the diquark is on its mass shell (see also 23|). 

As for the spin structure of the propagators, we have 
adopted the same scheme as in Refs.pH [5S] where we 
considered both the S = (scalar) and S = 1 (axial vec- 
tor) configurations for the diquark. This allows one to 
obtain distinct predictions for the u and d quarks. How- 
ever, we assume a similar form for the scalar and axial- 
vector couplings (scalar-like), and we distinguish their 
different contributions by varying their respective mass 
parameters in the calculations. This assumption is in 
line with previous estimates |24H26| where it was advo- 
cated that the full account of the axial-vector coupling 
does not sensibly improve the shape of parametrizations, 
while considerably increasing the algebraic complexity of 
the various structures (see e.g. dZ]). We define F as the 
scalar coupling at the proton-quark-diquark vertex 



^(P-Ml)2' 



f/s being a constant ^. The quark-proton helicity ampli- 
tudes are defined as 

^A',A';A,A= j d^k^4>l,^j^,{k',P')4>xAk,P), (19) 



^ The choice of this vertex function is motivated by phenomeno- 
logical reasons - it allows for an easier fit of the form factors — 
and also because it agrees with predictions within the Schwinger- 
Dyson formalism. See discussions in 1241 1251 and in |23| . 
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with 



'/'A,A(fc,P) = r(fc) 



u{k,X)U{P,A) 
fc^ — 



Finally, the denominators are evaluated with the di- 
quark mass on shell, 



defining the helicity structures at each soft vertex. One 



has: 



9AA = 9ax- 



for A = A, 



for A = -A. 



We list the separate structures appearing in Eqs.(|19l 



d^k^^ 


b\+{k',P')q 




(20a) 


d^ki^(, 


b\_{k',P')q 


^+~{k,P) 


(20b) 


d^k^c, 


t>*-+{k',P'h 


i++(fc,P) 


(20c) 


d^k^(, 


b*++{k',P'h 


i_+(fc,P). 


(20d) 



k'^-m^ = XM"^ - 



X 

l^X' 



1 - a: 



(21a) 



m2 = ^4a/2-^^m2 



1-C 1-a: 



X 



rn? — 



1-C 



1-C 

1 -X 
(21b) 



To extract the GPDs we calculate the convolutions in 



Eqs.(15l using the expressions for the g and A functions 



evaluated above. We obtain 



1-a: 
1-a: 



(m + MX) {m + + k_L • k_L 

{k^ - MlYik''^ ~ MlY 
(m + MX) (fci + ik2) - (m + Mff^) [ki + i/sa) 

(fc2-M2)2(/!c'2_M2)2 



l-C/2 ^(l-C/2)yW 



1 



yW(l - C/2) 2Af 



(22) 



E (23) 



1-a: 
1-a: 



[m + MX) (m + Afff^ 



-C2/4 



(fc2_Af2)2(fc'2-M2)2 l-C/2" ' (l-C/2)^/W 

{m + MX) {ki+ih) + (m + M^^ (fci +^^2) (^/2 



E 



Ai + iAj 



(fc2- A/2)2(A:'2_M2)2 



yW(l - c/2) 2A/ 



(24) 



E (25) 



from which the following forms for H, E, H, and E can be derived 



H ^ N 



(m + AfA:) ( m + ^^^^-^ ) + k_L • k_L 



(fc2- A/2)2(fc'2_M2)2 



4(W)^' 



(26) 



E = Af 



I- X 



-2M{1 - C) 



(to + A/X)^2 [m + M- ' 



1-cy Ai 



(fc2-A/2)2(fc'2„y\,j2)2 



(27) 



(m + AfX) ( m + M^^^—^ ) - k_L • k_L 



{k^ - Ml)^{k'^ - Ml)^ 



4(1-0 



E 



(28) 
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E = Af 



l-C/2 
I- X 



4M(1 - C) 
C 



(m + MX) 



k ■ A 
AT 



m + M- 



X-C\ k±-A 



(29) 



where is in GeV*. The integrations over cPk± yield 
finite values for the amplitudes in the limit A — > and 
C (see Appendix |A|. 



C. Reggeization 



It was noticed in Refs. [Ml [5S] that the low X behavior 
of the GPDs in the forward limit necessitates an extra 
factor of the type w X~°' in order to adequately fit cur- 
rent DIS data. We reiterate that, for the off- forward case 
this factor is important even at intermediate/large values 
of X and C because the dominant behavior of the GPDs 
at low X determines the nucleon form factors. In other 
words the lack of such a term hinders a good fit of the 
form factors. 

The Regge term indeed can be seen as originating from 
a generalization of the diquark picture in which the mass 
of the diquark, is not fixed, but has a spectral distri- 
bution. As first shown in Ref. [S^ in a simple covariant 
model for a pdf, choosing a spectral distribution of the 
form pii{M^) (X M^^^^ gives rise to the pdf behavior 



-a(0) 



upon integration over all Mj^. A more physical 



picture of the nucleon might have two contributions. 



p(Mi) 



(30) 




M„2 (GeVO 



FIG. 3: A plausible diquark mass distribution, the spectral 

1 

3 • 



density with Regge behavior M^"''''' with a(0) 



wherein pa is sharply peaked at a fixed value of the di- 
quark mass = Mx, and it yields the usual "fixed 
mass" diquark model. pR smoothly extends to large 
values and, with the Regge power behavior, upon inte- 
gration over , will yield a Regge type behavior largely 
dictated by the shape of the spectral function itself . 

We will illustrate this "reggeization" process by con- 
sidering the spin independent GPD H{X,0,0) = fi{X) 
as a function of a continuum of diquark masses. Aside 
from overall constant factors, the exact expression for the 
forward limit for the mass Mx, obtained from Eq.( 26 1, 
is given by 



fl"(X,0,0) cx 



2{mg + XMf 



[Ml + ^^(Af2 - XAf2)]3 
X 



[Ml 



X 



{Ml - a:m2)]2 



(31) 

6X3 



Multiplying this expression with pulMl) oc M^^^^ and 
then integrating over all diquark masses from to oo 
gives the analytic result 



dMlpR{Ml)H{X,0,0) ^ X-"(0)-i, 



(32) 



for X — > 0. In practice we would not integrate from 
zero mass and we would cut off the integral at the max- 
imum mass allowed by the kinematics. The fixed mass 
term pg{MV) would give back the unintegrated expres- 
sion Eq.(|32[) which goes to a constant or X^ as X ^ Q. 
A plausible form for the spectral density of the diquark 
is shown in Figure [3] 

Our main conclusion is that we can generate the Regge 
behavior and still be consistent with the diquark model. 

At this point this discussion does not consider the t- 
dependence, although it is plausible to incorporate the 
Regge trajectory form a{t) = ao + ct'^j while includ- 
ing the ^-dependence of the diquark model. Including 
the skewedness, ^ 7^ is more complicated, because 
of the distinction between the DGLAP and ERBL re- 



gions. While detailed calculations using Eq.(30) will be 



presented elsewhere, in the present analysis we adopt a 
factorized form of the Regge term which has similar fea- 
tures as the Reggeized diquark model, but is more in line 
with current parametric forms of parton distributions. 



(33) 



In this form only the term X " can be considered a 
proper Regge contribution. The term a'{X) = a'{l—XY 
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is constructed so as to guarantee that upon Fourier trans- 
formation in A , one obtains finite values for the nucleon 
radius as X — > 1 j^S]. This prescription is obtained for 
C = 0, or t = -A5_ [H]. In order to extend it to C 7^ 0, 
an additional term in the exponent, /3(C), is introduced. 
The physical motivation for this term is that it effectively 
accounts for the shift between the initial and final pro- 
ton's coordinates that occurs when Fourier transforming 
GPDs at C 7^ [Sn]. 

As we will show in Sections |III| and |IV[ two forms of 
/3 are suggested by the behavior of the available DVCS 
data, 



From this expression one defines 



/3/(c) = /^.^ _ ^ 

/3//(C) = /3C^ 



(34) 
(35) 



The effect of these terms is to allow for a data driven 
change in the slope in C of the GPDs, most likely an in- 
crease, with respect to the diquark model predictions. 
Since most DVCS data so far appear as asymmetries 
given by ratios of different cross sections combinations, 
it is difficult to determine precisely the C behavior of the 
CFFs and GPDs. It is therefore indispensable in future 
experimental analyses to provide absolute cross sections, 
as already done for the set of data provided by Hall A. In- 
troducing directly DVCS data to determine the behavior 
of our fit is an important step that distinguishes our anal- 
ysis from other ones in that it helps establishing the main 
trends of the multi-variable dependent data. The treat- 
ment of multi- variable dependence characterizes analyses 
aimed at extracting GPDs from data. What we suggest 
here is a bottom-up approach where experimental evi- 
dence is used to constrain the various theoretical aspects 
of the GPDs behavior, eventually giving rise to a com- 
plete model. 



D. Crossing Symmetries 

GPDs observe precise crossing symmetry relations. In 
order to discuss these symmetry properties we first intro- 
duce the so-called "symmetric system" of variables {x, ^}, 
where 



k+ + k'+ _ X-_C/2 
P+ + P'+ ^ 1 - C/2 



2A^ 



P'+ 2-C 



We also introduce the quark labels for all four chiral 
even GPDs, Fq = {Hg, Eg, H^, Eg}. By analogy with 
DIS, we define Fg{x,£^) in the interval —1 < a; < 1, with 
the following identification of anti-quarks, 

Fgix,0 = -Fq{x,0 x<0. (36) 



F^x,^), 



(37) 
(38) 



where Eg is identified with the flavor non singlet, valence 
quarks distributions, and J2q Fg with the flavor singlet, 
sea quarks distributions. Eg and E^ obey the symmetry 
relations 



Eg{X,i) = Eg-{-X,0 

Eg+ix,0 = -F+i-x,0- 



(39) 
(40) 



In DIS the commonly adopted Kuti-Weisskopf model en- 
sues [31] by which all distributions are evaluated at pos- 
itive X. 

In the off-forward case crossing symmetries are impor- 
tant for the evaluation of the CFFs deflned in Eq.([2|. 
The Wilson coefficient function, in fact, also obeys cross- 
ing symmetry relations. 



C^{x,0^ 



± 



1 



X — ^ + ie X + ^ — ie' 



(41) 



so that 



■Hg = f dxC+{x,OH+{x,tt) (42) 
Jo 

Uq ^ f dxC-{x,OH-{x,^,t), (43) 



similar relations hold for Eg and Eg. In the non- 
symmetric system of variables adopted throughout this 
paper, the axis of symmetry is shifted to X = C/2. More- 
over, X e [-1, 1] ^ X e [-1 + C, 1] , and x = ^ X = 

0, X = c ^ X = c. 

As explained in the Introduction, the validity of a sim- 
ple handbag based partonic interpretation of DVCS in 
the ERBL region has been recently questioned. The 
safest choice for a parametrization in the ERBL region is 
therefore to adopt a "minimal" model that accounts for 
crossing symmetry properties, continuity at the crossover 
points (X — and X — (^), and polynomiality. A possi- 
ble form is obtained by parametrizing the crossing sym- 
metric and anti-symmetric contributions, respectively as 
follows 

H^^^iXX) = a-(C)X2-a-(C)CX + ff(C,C)(44) 
H+^^{X,C) = a+X'^ -a+CX^ + cX + d (45) 

where a-(C) = 6(Ci?(C,C) - 2SEB.BL)/e, Serbl being 
the area subtended by Hg in the ERBL region, 



ERBL 



dXHg{XX,t) 

Jo 



dX (46) 
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FIG. 4: (color online) Crossing symmetric, ff"*", and antisym- 
metric, , contributions to the GPD aX C, = 0.18, and 
t — tmin- In this example, was taken equal to zero in the 
DGLAP region. 



Serbl appears in the definition of a~ multiplied by a 
factor of 2 because of the crossing symmetry pro perty 
for the areas subtended by H~ and Hq. In Eq.(45 1 a+ is 
a free parameter. This choice of parameters gives 



H = at X = the antiquark component in the 
DGLAP region being taken to be equal to zero in this 
phase of our analysis. Notice that Hq and Hq are not 
required to obey crossing symmetries. They are obtained 
by construction from Eqs.( 44p5 ). An example describing 
the symmetric and antisymmetric components of -ff' is 
given in Fig|4] 



III. RECURSIVE FIT: NUMERICAL 
EVALUATION OF GPD PARAMETERS FROM 
INCLUSIVE MEASUREMENTS CONSTRAINTS. 



Here we describe a recursive fitting procedure to ex- 
tract the chiral-even GPDs from available DVCS data. 
Our fit uses the parametric forms (we omit the quark 
labels for simplicity) 

FiXX,t)=MGfj-jX,Ct)R;''^'iX,Ct) (47) 



and i?p was given in Section II C 

It should be remarked that our new parametric form 
follows from the one used in Ref . [24l |25] , while present- 
ing several important differences. We have first of all, 
completed a thorough analysis of the spin components 
of the various GPDs, both in the unpolarized and polar- 
ized sectors, thus releasing the assumption of a simplified 
quark-proton vertex structure made in [24, 25j, and ex- 
tending our analysis to the much needed H and E func- 
tions. The more careful spin treatment also results in a 
different shape for E which in 24, 25 closely followed H 
by construction. 

The most important features of our new parametriza- 
tion are summarized below: 

i) we consider only configurations for a spin 1/2 quark 
and a spin diquark. The flexibility in shape contributed 
by considering a spin 1 diquark was in fact not sufficient 
to allow us to model e.g. the rise at low X. We therefore 
opted for keeping the Regge term as in [24l [25]. This 
can in fact be derived from a "reggeized" version of the 
model as explained in Section [Til 

ii) we model all chiral-even GPDs, and we present for the 
first time parametric forms for H and E, besides new 
evaluations for H and E. Our analysis applies to the 
intermediate xsj, multi-GeV regime which is domi- 
nated by valence quarks in the DGLAP region. Only u 
and d quark flavors are considered. 

in) we perform a recursive fit in which parameters are 
evaluated orderly, from imposing constraints from DIS 
experimental results first, then from the elastic form fac- 
tors, and eventually including DVCS data directly. This 
procedure affords us a better control on: i) the number 
of parameters that are necessary to constrain the GPD 
multi-variable problem; ii) the fit's variants as new data 
are inserted. 

All parameters obtained from the DIS and elastic con- 
straints are given in Table |T] They correspond to the first 
two steps of our fitting procedure. While we address in 
detail the impact of the GPDs H, E, H on the description 
of available DVCS data, we postpone the discussion of E 
to a dedicated analysis in [32j. As we explain in what 
follows E contributes to DVCS observables multiplied by 
a factor x bj , or ^ « xbj /{2 — XBj). As experimental data 
on DVCS target asymmetries accumulate, it is important 
to clarify that the pion pole contribution to E scales as 
1 /x Bj only in specific models like the chiral soliton based 
factorized form described in Ref. [33]. While this factor- 
ized form is a convenient model, the ^ singularity is not 
required by the general analytic structure of the GPD. In 
fact, in our evaluation E is estimated to be suppressed 
by a factor < 0.1 at HERMES kinematics. 

The first set of experimental constraints is given by 



where F = H,E,H,E; the functions G^A ^ are the 



covariant diquark contributions from Eqs.( 26|27|28|29 |, 



^ In Appendix|B]we present additional parametric forms that are 
more practical for applications and numerical calculations. 
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the valence contribution to the inclusive DIS structure 
functions, 



i7«(X,0,0,Q2) 
i?«(X,0,0,Q2) 



/f(X,Q2) = g,(X) (48a) 
gHX,Q^) = Aq,{X) (48b) 



representing the forward limit of Eq.(47l (we have re- 
stored both the quark labels q — u,d, and for the 
dependence), /f and gf are obtained from DIS data. We 
do not use directly experimental data in this phase of the 
analysis, but we perform a fit of the valence components 
of existing parametrizations. The fit was performed sim- 
ilarly to Refs.[2l[2g. By inspecting Eqs.Q, ([28]) and 
(471 one can see that for t = and C — 0; the only pa- 



and 



rameters that enter are: Mx, M\, m m 

We fit separately the unpolarized, H, and polarized, H 



GPDs. For H, an additional parameter, is fixed by the 
baryon number sum rules constraints, dXuy{X) = 2, 

and Jp dXdv{X) = 1. Therefore in our first step we have 
four parameters per distribution, per quark flavor giv- 
ing a total number of parameters consistent with what 
obtained in recent PDF parametrizations e.g. |34H36j . 
As already noticed in the diquark model based 

parametrization corresponds to a low initial scale, Q^. 
Parametric forms are then evolved to the of the data 
using LO Perturbative QCD (PQCD) evolution equa- 
tions [37] . Additional parameters not shown in the Table 
are the initial value of the perturbative evolution scale, 
= 0.0936 GeV'^, and the parameters /3 appearing in 
Eq.(|34]), /? = 10, and Eq.(|35]), /3 = 1.5. These were fixed 
by implementing directly DVCS data in our fit, as we 
will show in Section IIVI 



Parameters 


H 


E 


H 


E 


m„ (GeV) 


0.367 


0.367 


2.479 


2.479 


(GeV) 


0.583 


0.583 


0.467 


0.467 


MX (GeV) 


0.963 


0.963 


0.909 


0.909 


Ou 


0.222 


0.222 


0.218 


0.218 




2.443 ± 0.063 


4.582 ± 0.128 


1.758 ± 0.839 


5.549 ± 0.519 


Pu 


0.6649 ± 0.0268 


1.465 ± 0.031 


0.558 ± 0.468 


0.420 ± 0.069 


JVu 


1.468 


1.468 


0.0343 ± 0.0033 


4.882 ± 0.636 


md (GeV) 


0.0850 


0.0850 


1.211 


1.211 


M| (GeV) 


0.841 


0.841 


0.699 


0.699 


Mi (GeV) 


0.7592 


0.7592 


0.836 


0.836 




0.0378 


0.0378 


0.0417 


0.0417 




1.777 ± 0.021 


0.0516 ± 0.0026 


1.489 ± 0.629 


4.791 ± 0.316 


Pd 


0.114 ± 0.015 


-10.147± 0.681 


1.032 ± 0.552 


0.248 ± 0.060 




1.023 


-2.368 ± 0.160 


-0.0768 ± 0.0068 


-17.414 ± 1.815 



TABLE I: Parameters obtained from om recursive fitting procedure applied to Hq, Eg, and Hq, q = u,d. 
and Oq were obtained in a first phase by fitting valence quarks PDFs from DIS experimental data, a'q and pq were obtained 



Ml, 



Ml 



subsequently, by fitting the proton and neutron elastic form factors 
of the additional parameter /3 in Eq.(33l is /3 = 10. 



and the axial form factors 



In the fit for H we use a similar scheme as in current fits 
[311 [31], where 



XAq,{X, Q2 



^NqX- 



^'•Xq,iX,Q^). 



We left the mass parameters Mx and M\ fixed as for the 
unpolarized case, while we varied a and m. By letting 
the latter vary, we obtain the effect of the extra term 
oc (1 -t- 7^) introduced in [S^ for the LO fit. Figure 
[5] shows our curves for Hu,d{X, 0,0) at the rather high 
value of = 25 GeV^, in order to test the stability 
of our fit with PQCD evolution. Other available PDF 



parametrizations from quantitative fits are also shown in 
the figure. 



In the second phase of our fit we impose an additional 
set of independent experimental constraints from the nor- 
malizations of the chiral even GPDs to the nucleon form 
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FIG. 5: (color online) GPDs H^iX, 0, 0) (top) and HaiX, 0, 0) 
(bottom), evaluated at = 25 GeV^, compared with current 
LO parametrizations [34H36] . 



factors, 







(49a) 






(49b) 






(49c) 






(49d) 



where -F' (t) and ^2(0 the Dirac and Pauh form fac- 
tors for the quark q components in the nucleon. G\{t) 
and G'^p{t) are the axial and pseudoscalar form factors. 

Notice that the GPD E is constrained by the pseu- 
doscalar form factor of the nucleon through Eq.(49d|. 



When the covariant or light front diquark spectator 
model is applied to calculating i?(X, C,t), there is no 

kinematical singularity. The combination y- + 

A ^__|_+ vanishes as C — for any X and t. This ap- 
pears as a restriction on the GPD in the DGLAP region, 
X > C,. Requiring polynomiality leads to the sum rule in 
the ERBL region 



dX 
l-C/2 



E{XX,t)^Gp{t)- 



dX 



C l-C/2 



E{X,Q,t) 
(50) 



Since the integral in the DGLAP region is finite for any 
C and does not have a pole at ^ = 0, the ERBL region 
integral will not either. This will be true of any spec- 
tator model wherein there are no kinematic singularities 
introduced. In the diquark spectator approach that we 
use, the t dependence of the pion pole in the form fac- 
tor can be reproduced while satisfying the sum rule in 
) by a suitable choice of mass and 
This corresponds to a dual picture 



Eq. |50| (for small \t\ 
parameters 



"Regge" 

- a ^-channel pion pole emerging from an integral over an 
s-channcl diquark pole. In summary, we reiterate that 
the GPD E enters the target asymmetry always multi- 
plied by C,{xBj)^ so that it contributes only weakly in the 
HERMES kinematical region. 

Isospin decompositions allow one to relate the quark 



the 



form factors to experimental measurements of -^{(2) 
Dirac (Pauli) form factors for the proton (neutron), re- 
spectively. Ga and are the isovector and isoscalar 
components of the axial nucleon form factor, and Gp. 
We used the same selection of data as in Ref.pij for 
the nucleon electromagnetic form factors (see references 
in |40| ) . The resulting parameters are given in Table 
|lj More recent data |1T] are now available that show a 
milder slope of the electric to magnetic proton form fac- 
tors ratio at large However, these do not largely affect 
our fits that are limited to the — t << w 2 — 3 GeV^ 
region. Ga is obtained from the global average of neu- 
tron beta decay, and neutrino scattering experiments (see 
Ref. [42] and references therein), 



GA{t) 



9A 



Ml 



(51) 



with gA = 1.2695±0.0029, and Ma = 1.026±0.021 GeV. 
Gp is notoriously dominated by a pion pole contribution, 
a small non-pion pole component being also present. We 
used the experimental values displayed e.g. in Ref. |43j . 
A more thorough discussion of this form factor will be 
given in [521 . 

As a result, for each quark flavor and GPD type, using 
the constraints above, we can determine the additional 
parameters, a', p in Eq.(47l and the normalizations M 
(Eq.(47l and Table [l|. The number of parameters used is 
consistent with the one used in fits of the nucleon form 
factors data. These require four parameters for 
G'li, respectively, and two for Gg. [iPllil] . 

The GPDs H^JX,0,0), 0, 0), i7„^<i(X, 0, 0) 

are shown in Figls] both at the initial scale, Q^, and 
evolved to = 2 GeV^. A comparison with results on 
H'^ and E'^ from [24l [25] at the scale Q^, is also shown. 
In Fig[7]we show ff„_d(3f, C, t; Q^) evaluated at = 2 
GeV^ and for a variety of ranges in C = XBj and Q^. 
In Fig[8] we show the working of the property of poly- 
nomiality. This is, in a nutshell, a direct consequence 
of extending the Operator Product Expansion (OPE) to 
the off-forward case [5|[1S], according to which the Mellin 
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FIG. 6: (color online) GPDs i^,(X,0,0) = {Hg,Eq,Hq}, for 
q — u (left) and q = d (right), evaluated at the initial scale, 
Qo = 0.0936 GeV^ and at = 2 GeV^ respectively. The 
dashed lines were calculated using the model in Refs.|24l 125) 
at the initial scale. 



FIG. 7: Hu{X,C,t;Q'^) evaluated at = 2 GeV^ Each 
panel shows Hu plotted vs. X at different values of ^ = 
0.18,0.25,0.36,0.45. For each value of ( several curves are 
shown that correspond to a range of values in ~t from tmi„ = 
-M^CVCl - C) to 1 GeV^ 



moments of GPDs read (see also [25]) 

= A,,o{t) + A..2(<)(20' + - + A„,„_i(t)(20"-(52) 

where n > 1, and only even powers of ^ are included. 
Similar results hold for Eq,Hq,Eq 46 ; for n > 1 the 
equation is also dependent. To illustrate polynomial- 
ity, the moments of were plotted vs. ^ at the initial 
scale, Qo, for different values of t, and for n < 5 (the 
trend shown in the figure holds for even larger values of 
n). The two sets of curves represent the calculation using 
the parameterization from this paper, and the results of 
a polynomial fit in . Clearly, our parameterization sat- 
isfies the property of polynomiality although this cannot 
be inferred directly from the functional form in Eq. ( 47 1 . 



To summarize, we proceeded through two steps of our 
recursive fit. We used the flexible model described in Sec- 
tion |ll] to first fit the nucleon unpolarized and polarized 
PDFs, respectively, by using all parameters that enter the 
expressions a,t t — and C = 0. We subsequently fitted 
the nucleon electromagnetic and electroweak form fac- 
tors using the additional parameters that enter at < 7^ 
and C = 0- At this stage of the analysis we established 
that in order to obtain GPDs that are constrained by a 
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FIG. 8: The polynomiality property in our parametrization. 
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set of parameters which is consistent with the ones used 
for independent fits of the DIS structure functions, and 
of the nucleon form factors, a total number of 4 (DIS) 
plus 3 (Elastic) parameters per quark flavor, per GPD 
is needed. The number of parameters is, in practice, re- 
duced because of the physically motivated ansatze intro- 
duced in our approach, as one can see by inspecting the 
values in Table [T] Physical assumptions are both useful to 
understand the trend of data and at the same time they 
are known to introduce a bias. This aspect is well known 
to affect all hard processes multi- variable analyses, and 
it can be perhaps overcome in suitable neural network 
based approaches [371 EH] ■ 

The third set of experimental constraints is given by 
DVCS-type data. In our analysis we use strictly DVCS 
data that are cleaner from the theoretical point of view, 
as compared to Deeply Virtual Meson Production. By 
fitting the CFFs that are functions of C, t and Q^, we will 
be able to provide a constraint on the otherwise elusive C 
dependence of the GPDs. To understand how many extra 
parameters are needed in extending the fit to its third 
phase, we devised a procedure explained in the following 
Section. 



Sections ITT] and III In order to proceed we first introduce 



the helicity formalism. This method allows us to obtain 
a physical interpretation of the various structures, and 
azymuthal angular dependences involved in terms of the 
photons helicity states. 



A. Observables 

The observables included in our fit are from all the 
DVCS measurements that were available to us, to date. 
These are: the cross section for unpolarized electron 
scattering, da/d^ [321 120] j the Beam Spin Asymme- 
try (BSA), Alu [311 ED, the Beam Charge Asymmetry 
(BCA), Ac [52], and the Transverse Spin Asymmetries 



(TSAs), A^^<^^, and A{,j, [5l [S3]. The cross section 
for scattering of an electron/positron beam with polar- 
ization /i off a proton target is evaluated considering the 
sum of the DVCS and Bethe-Heitler (BH) amplitudes. 



d^a^- 



dQ'^dx Bj dtdipdcl) 



r I r, 



BH 



DVCS 



N 



(53) 



where the factor P is given by 



IV. IMPLEMENTATION OF DVCS DATA 

We now discuss our procedure to extract GPDs from 
available DVCS data. Fully quantitative studies were 
performed in [161 117j where a number of observables were 
fitted, that were obtained from experimental data for the 
process ep — ep^ from both Hermes and Jefferson Lab. 
In order to have a consistent comparison, we included 
in our analysis a similar set of observables. Below we 
list their expressions in terms of the CFFs displayed in 



XBjV 



16^2 g2(2AfxBjei) VT+^ 

with y = v/ei, ei being the initial electron energy, and 
V = 61 — 62 the momentum transfer; 6 = AM^x\^/Q'^] (j) is 
the (azymuthal) angle between the hadronic and leptonic 
planes in the frame where the virtual photon's momen- 
tum is along the z-axis [10]. The unobserved helicities 
have been summed over, implicitly. ^ The various ob- 
servables that we consider are written as 
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da 



A 



LU 



Ac = 



iDVCS 



A' 



d^CJ^ d^CT^ 



d^ 




d^a"^ 






d$ 


d*<j+ 













= r [I Tbh P + I Ti3vcS I' + (/^ + /^)] , 

r(/i- ~ /^) 

2da/d^ ' 

r(/^ + /^) 

2da/d^ ' 



5, 



r I Tj^/gg p 

2da/d^ ' 



TP 



5, 



2da/d^ 



(54) 
(55) 

(56) 
(57) 
(58) 



where d$ = d(f>dxBjdtdQ^; the superscripts +(— ) refer to 
the beams' charge, f (i) are for oppositely polarized elec- 
tron beams, the subscripts (=>) represent the trans- 
verse target polarizations, corresponding to the angles (j)g 
and (ps + TT, respectively. The subscript TP follows the 
notation of (TU] for transverse polarized target; the sub- 
script LU is for a longitudinally polarized beam, L, and 
an unpolarized target, U, while UT is for an unpolarized 
beam, and a transversely polarized target. 

An expression for Tbh, the amplitude for the Bethe- 
Heitler (BH) process is given in Ref.|10]. Here we write 
the amplitude in helicity basis in order to facilitate the 
expansion of the observables in bilinear products of am- 
plitudes and GPDs. The basic form of the BH amplitude 



r; 



,/i,A.A',A' 



BH 



2 "^A.A' ' 



where 



and 



"^A,A' 



C7(p',A') [i^i(A2)7'^ 



2M 



U{p,K). 



(59) 



(60) 



(61) 



The DVCS amplitude for scattering of a lepton with spin 
h is given by 



pAA'.A; 

^DVCS 



1 



= u{k2,h)'yf,u{ki,h)—Tj^''^,e,y 



* A' 



(62) 

where it(fci(2),^) are the initial and final lepton spinors; 
the hadronic tensor, Tj^\,, was defined in Sectionjllj and 



£[/ ^ is the outgoing photon polarization vector. Tdvcs 
can be expressed in terms of helicity amplitudes by con- 
sidering the following expansion on the polarization vec- 
tors basis EH ESI , 



T 



AA'A' 



DVCS 



A+/+.A;A1,,A 



A /-,A;A;,A' 

^3 /o,A;A' ,A' 



(63) 



with 



A+ = 




T^u{k2,h)[^i±i-f2\u{ki,h) (64a) 



A. 



2e 

Q2\/r^ 

u{k2,h)^3u{ki,h) 



(64b) 



The / amplitudes were given in terms of CFFs in Sec- 
tion II. At LO, by disregarding the longitudinal photon 
polarization, the only amplitudes that were found to con- 
tribute are: /++,++,/-+-+,/++,+-,/-+,— ■ As a con- 
sequence, the only term contributing to the unpolarized 
term, | Tdvcs corresponds to the transverse cross 
section, dar/dt, and it is given by 



T 



1 



1 



DVCS 



Q2 2(1 - e) 

I /++,+- P 



/++,++ r + 1 /-+,-+ 

I /-+,- n 



(65) 



In terms of CFFs, 
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1 



1 



^DVCS — 



2(1-6) {2-XB,] 

/ ^2 



t 



\-X 



Bj 



4M2 



4(1 - XBj)in*n + n*n) - xl^irn + %*£ + rn + %*£)- 

I 



Eq.(66l is analogous to term Cq^^^ in the expansion 
given in [TD] . Note that the sum over all the unobserved 
helicities is implicit. 

An analogous decomposition into the virtual photon 
polarization basis for the Bethe-Heitler amplitude of 
Eq. ( 59 I has the form 



/i,a,a',a; _ + 



^/i.a;' 

The hadronic amplitudes are 

'^il.AiA'^A 



(67) 



J± 



1,A:A' = -A 



= A 



2M^\ - ^2 



0A,± 



{F^{t)+F2{t)), (68) 



(66) 



where terms of order (Aj^/p+)2 and {Ad/p'^Y were 
dropped. 

The lepton tensor for the Bethe-Heitler amplitude can 
be calculated from Eq. ( 60 1 using the relation 



^i^YY = g'"'!" + g""!^ - ^^7" - le^'^'T^Ts. (69) 



Momentum conservation gives the exchanged virtual 
photon momentum as = k2 + q' — ki = q' — q = 
q' ~ {p — p') ~ q' — A. It can be seen that the Dirac 
algebra is reduced to evaluating a single 7*^ or 7'^7^. 

The interference term is given by, 



/ = 



E 



h,AA'A' /i.AA'A', 



BH 



DVCS 



' h,A.A',A' /i,AA'A' 



DVCS 



BH 



(70) 



h,AA'A' 



We consider the following expansion over helicity states of Tg^ 
lepton helicity h, 



/!.,AA'A-, 



and T, 



/i,A,A',A' . 



DVCS 



in terms of / 



A,,A:A' .A 



/ for fixed 



( 



-B^_l_*i^-l; + ,+ + Jo: + ,+ + iL*! J+l; + ,+ ^h^) /+. 



+;+,+ 



where the various helicity dependent terms are defined in Eq. ( |60|68|63| . 

Eq.(70) then has the following structure analogous to the leading terms in the expansion in '101, 

/ oc Cg + c( cos (j) + hs{ sin 



(71) 



(72) 
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with coefficients given by 



, (2-2?y)= 
1-2/ 



2 - XBj 4M2 



(73) 



c[ = -8{2-2y + y^)K 



Fr&eH 



XBj 



2 - XBj 



{Fi + F2)^eH 



4M2 



F2^e£ 



(74) 



82/(2 - y)K 



Fi^mH + 



2 - XBj 



{Fi + F2)^m'H 



4M2 



F2^m£ 



(75) 



The kinematical factor K is, at leading order in —t/Q'^ 



K = 



-{l-XBj)(l-y-^ 



(l + e) 



1/2 



1/2 



with train = -X%j/{1 - XBj)M^. 



For a practical approach wc streamline both the kine- 
matical dependence and the GPD content of the various 



observables. By keeping the leading terms in —t/Q'^, and 
in XBj, we obtain 



d(Te- 
./'I' 

Alu = 
Ac = 



aDVCS 



r (B((/>) + C(0) COS0) 

A((/)) sin^ 
B{<j)) + C{(l)) coscj) 

D{(j)) cos(j) 
B{(t>)+C{(j)) cos (I) 
Ejcj)) sin(0 - (t>s) 
B{(j>) + C{(j)) cos(/) 
{F{^) + G{(f)) cose 



■) sin((/) - 4>s) + F[{4>) cos{(j) ■ 
B{(t)) + C{^) COS0 



(76a) 
(76b) 
(76c) 
(76d) 
(76e) 



where, 
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A 

B 

C 

D 
E 
F 

G 

H 



Ki{(j)) Fi ^mn 



XBj 



-{Fi + F2)^mn 



2 - XBj 

Ci"i(t)) + K%cs {^e^'H + ^m'^n + ^e^H + + Relief)) ( Fi^eH + 



Ki{cj)) Fi ^eU 



2 - XBj 

^^^^{F.+F^neu) 

2 - XBj J 



{Fi+F2Wen 



= K 



UT,0 
DVCS 

I 



KY^'°i^) (^{2~XBj)F,^m£ 



= K 



K 



I VP) 

UT, 1 



V4Af2 



— (2 - XBj)Fi 3m£ - 



t 1 



XBj 



4M2 



X Bj Fi 5m£ 



]VP2-XBj 

t 1 — XBj 



F2'^mH 



M2 2 - XBj 



Fo '^roH 



XBj 



XBj 



(77a) 

-{Fi^F2Wen]{71h) 

(77c) 

(77d) 
(77e) 
(77f) 

(77g) 
(77h) 



r 



C^^((/)) enter the BH cross section [10 . The fac- 
tors K are kinematical coefficients which depend on 
t,XBj,y,Q^,^', for the BH and interference terms, they 
depend also on </> due to the BH propagators. 



B. Fit Results 



We present resuhs from our fit including all parameters 
that were fixed using the reggeized diquark parametriza- 
tion described in Sections |TTj |HI[ and displayed in Ta- 
ble [ij We introduce additional flexibility through extra 
parameters entering Eq. ( 33 1 in order to constrain the C 
dependence from all available DVCS data. Although in 



principle as many (3 parameters as the number of flavors 
and GPDs could be introduced, given the small data set 
presently available, we use one value of (3 for all GPDs. 
More accurate studies including flavor and GPD depen- 
dent /3 parameters will be conducted as more data be- 
come available. In Fig|9]we show the real and imaginary 
parts of the CFFs, 'H{C,t), appearing in Eqs.(77|. Sim- 
ilar results are obtained for E, H, and E. One can see 
that the slope in ( flattens out a.s —t increases. 

Our fit uses the two currently available sets of data, 
from both Hall A and Hall B collaborations at Jefferson 
Lab, and the Hermes collaboration, respectively. Since 
the kinematical ranges covered by the two experiments 
only partially overlap, in this first step we start by fitting 
the Jefferson Lab set, and subsequently extend the results 
of our fit to predict the behavior of the Hermes set. 
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FIG. 9: (color online) Real and imaginary parts of the CFFs, 
HiCt), Eqs.([77|. The CFFs are plotted vs. xsj = C, for 

" = 2 GeV^ 



different values of i, at 
obtained for E and H. 



Similar results are 



1. Jefferson Lab Data 

The DVCS data on Alu from HaU B [5l|, and on the 
sum and the difference of the beam polarization cross 
sections from Hall A |49j were implemented directly in 
our fit. This allows us to constrain the parameter /?. In 
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FIG. 10: Beam spin asymmetry, Alu{4> = 90°) in 12 of the 
XBj and bins measured in Hall B [51 . The data points 
were extracted by fitting j4lc/(0), however they do not repre- 
sent the uncertainties reported in the experimental analysis. 
The second panel from the top includes also data from Hall A 
|49| . The full circles, open circles, and triangles represent data 
in similar xsj and t bins, but at values slightly displaced 
around the value reported in the legend. All curves were ob- 
tained at the kinematics displayed in the figure. Dashed lines: 
results from the fit using only the PDF and form factors con- 
straints as from Table |l] The full lines rep resent the effect of 
introducing the ^ dependent term, Eq.(33l, in the numerator 
of the asymmetry only. 



Figs[TO] and [TT] we show the results of our fits using the 
quantity Am{90°) = A/B, Eq.(|76b[). The experimental 



points in the figure were obtained by fitting in 
the 12 bins displayed in Ref.[5T] (figure 4). The statisti- 
cal and systematic errors were added in quadrature, no 
error correlations were considered (this is giving rise in 
our case to larger error bands, although the central points 
coincide with the ones in Ref.|51|). A similar procedure 
was used for the Hall A data that are also displayed in 
the figures. The dashed lines in both figures are a pre- 
diction of the fit using only the PDFs and form factors 
constraints. Clearly by taking only these constraints, the 
C slope of the CFFs is unconstrained and evidently off the 
data trend, as it can be seen in the larger ^ = xsj bins. 
The additional term in Eq.(33) can regulate this behav- 
ior. Two possible ways of implementing it are shown 
respectively by the full lines in the two figures. In Fig[TO| 
a multiplicative term was considered only in the numer- 
ator of the asymmetry, given by A. In Fig 11 the zeta 



FIG. 11: Experimental data: same as Fig ]10| Dashed lines: 
results from the fit using only the PDF and form factors con- 
straints as from Table |l] The effect of the dependent terms 
given in Eq.(34l, and Eq.(|35|l, is given by the full lines and 



the dot-dashed 



lines, respectively. 



dependent factor was introduced in the GPDs, and used 
to calculate the CFFs in both the numerator and denom- 
inator of the asymmetry. The effect of introducing such 
term in the GPDs gives a different dependence that can 
be ascribed to modifying both the real and imaginary 
parts of the CFFs. We show results using two different 
expressions for /3, Eq.(l34|) and Eq.(35l, in order to il- 



lustrate some of the subtleties that are involved in the 
extraction of the CFFs, and GPDs from the data. While 
the two expressions give almost identical results for the 
asymmetry, they impact the various terms. A, B, C in 
Eqs.(77) in different ways. With more data in hand, in- 



cluding a separation of the absolute cross sections, one 
will be able to perform precise fits of the behavior in C. 

In Figs 12 13 14|15 we analyze the effect of the differ- 
ent GPD components on the fit to Hall B data. Fig[T2| 
shows the contribution of the BH term and of the co- 



efficients A, and B from Eqs.(77l with the C dependent 



term from Eq.(34) (bullets), and without it (full curves) 



Fig |13| shows the separate contributions of the numerator 
and denominator, A, and B, Eqs.(77), to the ratio defin- 
ing the asymmetry ^1,^/(90°). B is given by the sum of 
the three terms displayed in the figure. Notice that the 
real part of the CFFs enters only T^yc's, which is a rel- 
atively small contribution. To understand its impact on 



the various observables, in Fig 14 we plotted results both 
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FIG. 12: Contribution of the BH term (dot-dashed hne) and 
of the coefficients A, and B from Eqs.(77l with the ^ depen- 
dent term from Eq.(34l (bullets), and without it (full curves), 
at = 1.7 GeV^ and xbj = 0.25. Similar results are ob- 
tained in other kinematics. 



FIG. 13: Separate contributions of the numerator and denom- 
inator, A, and B, Eqs.(77l, to the Alu{90°)- B is given by 



the BH term (dot-dashes) plus the sum of the [Tdvcs]'^ and 
interference terms (dashes). = 1.7 GeV^ and xsj ~ 0.25; 
similar results are obtained in other kinematics. 



2. Hermes data 



including and excluding the real part. 
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we show the effect of the GPD H 



Finally, in Fig 
on the fit. We confirm the result also quantitatively re- 
ported in Ref.[17] that DVCS data from an unpolarized 
proton target at Jlab kinematics are dominated by the 
contribution of the GPD H. The dashed curves in the 
figure were obtained by disregarding the contribution of 
H. E and E have also a very little impact on the data 
fit. 

In FigjlHjwe show the results of our fit vs. Hall A data 
These are given as the "sum" and "difference" of 
the two polarizations for the electron beam. Together 
with the data we also plot the results of a fit performed 
in [35] (yellow bands). All theoretical curves are shown 
with and without the C dependent correcting factor from 
Eq. ( [33| . Moreover, while we confirm that the sum, or ab- 
solute cross section part is dominated by the BH term, 
we also point out the importance of the contribution from 
the pure DVCS scattering, or | Tjjvcs at leading or- 
der in Q^. Also shown are the theoretical error bands 
for the asymmetry. In Fig |17| we compare our results to 
the experimental extraction of the imaginary part of the 
interference term coefficient in A^jj Eqs.(76b) and (77 1. 



The role of the CFF for H proves fundamental in deter- 
mining the slope vs. —t of this term. In order to illustrate 
this, we show a comparison with a previous calculation 
of the interference term where H was not included 1251 . 



In the second phase of our analysis, we use our fit re- 
sults to Jlab data to predict the quantities, Aljjj -^C, 
and AuT extracted at Hermes [52l [53] . Hermes data are 
provided as "coefficients" of the azymuthal angles de- 
pendent terms. The dependence of these coefficients on 
the various kinematical variables is sensitive to the set 
of approximations that one uses in the extraction, thus 
affecting quantitative analyses. In order to facilitate the 
comparison, and to once more show some of the sub- 



tleties involved, in Fig 18 we show the results of our fit 
for AlU vs. ~t, Q"^ , and xsj, respectively, calculated 
at each kinematical bin provided by Hermes [56) (curve 
denoted as "Hermes kinematics"), and at the nominal 
average values presented in each panel. It is interesting 
to notice that due to the correlation between x bj and 
in the data, different features arise when using the aver- 
age bin values. In the figure (lower panels) we also show 
the effect of disregarding the DVCS term in the denom- 
inator. Similarly to the Jlab results, the GPD that the 
data are largely sensitive to is H, the role of the other 
GPDs being marginal. As a concluding remark on Am, 
we notice that Jefferson Lab Hall B data seem to suggest 
a decrease of Ai^u with xbj- The curves in Fig 18 used 



Jlab data in the fit, and therefore they show a definite 
slope in xsj- 

In Fig 19 we show the different cross section compo- 
nents contributing to Alu(90°), given by A, and by the 
contributions from BH, Tpi/r:.?; and the BH/DVCS in- 
terference terms in B, Eqs(77). The curves in the figure 
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FIG. 14: Impact of difFerent GPD components on data fit. 
The effect of including (full curves) and excluding (dashed 
curves) the real part of the CFFs is shown for each component. 

— 1.3 GeV^ and XBj = 0.12; similar results are obtained 
in other kinematics. 
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FIG. 15: Impact of different GPD components on data fit. 
The effect of including (full curves), and excluding (dashed 
curves) H in the CFFs evaluation is shown. The contribution 
of E is always negligible; = 1.3 GeV^ and XBj = 0.12. 
Similar results are found in other kinematical bins. 



were calculated for the same kinematical bins as in Fig[T8| 
[56j . In order to discern the role of the various compo- 
nents involving GPDs from the kinematics, we also show 
"dry runs" of our code in the lower panels obtained by 
setting all GPD factors equal to one. The various kine- 
matical coefficients used in this analysis are written in 
Appendix [D| 



Finally, in Fig 20 and Fig 21 our predictions for Ac 
and are shown vs. —t. The agreement with Her- 

mes data is remarkably good within our theoretical error, 
despite we did not implement directly such data in the 
fit. 



V. CONCLUSIONS AND OUTLOOK 

In this paper we have presented a parametrization of 
the chiral even GPDs that is inspired by a physically 
motivated picture of the nucleon as a quark-diquark sys- 
tem with Regge behavior. The spin structure of each 
of the four GPDs is determined via the covariant quark- 
nucleon scattering amplitude, with a diquark exchange. 
The masses, couplings and Regge power behavior that 
set the scale for the dependence on the kinematic vari- 
ables, X,(^,t,Q'^, are determined via a recursive proce- 
dure. We fit the parton distribution functions fi and gi 
for the u and d quarks with H{X, 0, 0) and H{X, 0, 0) at 
a low scale. The electromagnetic form factors, Fi{t) and 
F2{t) constrain the first X moments of H{XX,t) and 
E{X,(,t). These first moments are constrained to sat- 
isfy polynomiality, thereby removing the C dependence 



0.1 
0.08 
0.06 
0.04 
0.02 
-0 
-0.02 
-0.04 

0.03 
0.02 
0.01 


-0.01 
-0.02 







DVCS 










Interference 






C^--.. Total 














Sum ((ib/GeV^) 













Difference (^b/GeV^) 







20 40 60 80 100 120 140 160 180 



FIG. 16: (color onhne) HallA data [49] for the "sum" (up- 
per panel) and "difference" (power panel) of the two electron 
beam polarizations. Shown are curves including the contribu- 
tion of the dependent factor from Eq.(33l (full lines), and 
neglecting it (dashed lines). All terms (DVCS, Interference 
and Total are shown for the "sum" graph. The yellow bands 
in both panels represent the error of the data fit. The green 



band in the asymmetry graph is the theoretical error from our 
parameterization. 
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FIG. 17: Coefficient C, Eqs.(72) extracted from Hall A data 
[491 . Shown are the contributions from the GPDs, H, E and 
H. 



All curves include the term in Eq. ( 33 
a previous prediction based on a simp 
and including H only [25] is also shown. 
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FIG. 18: (color online) Calculations at Hermes kinematics 
[521 |53l [56] . Shown is Alu{^^°) vs. —t, Q^, and xbj, respec- 
tively, calculated at each kinematical bin provided by Hermes 
|56| (curve denoted as "Hermes kinematics"), and at the nom- 
inal average values presented in each panel. It is interesting 
to notice that due to the correlation between xsj and in 
the data, different features arise when using the average bin 
values. In the lower panels we also show the effect of disre- 
garding the DVCS term in the denominator (dashed curves). 



FIG. 19: (color online) Cross section components contribut- 
ing to Alu{90°)'- a, BH (dot-dashed lines), the sum of the 
Tbvcs S'lid the BH/DVCS interference terms (dashes), and 
B, Eqs.(77|. The curves in the figure were calculated for the 
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same kinematical bins as in Fig 18 
the role of the various components involving GPDs from the 
kinematics, we also show "dry runs" of our code in the lower 
panels obtained by setting all GPD factors equal to one. 



and leaving only the t dependence. This requires fixing 
the parameterization of the ERBL region, X < C,, for all 
C, so as to satisfy a sum rule for the form factor. The same 
approach is used for the axial vector form factor and H . 
Similarly the pseudoscalar form factor constrains E, al- 
though we have not used that here (the contribution that 
is not dominated by the tt pole is not well known and is 
the subject of a forthcoming paper [32] )• 

In a previous paper Goldstein and Liuti [IH] have 
shown that the simple parton interpretation of the ERBL 
region is dubious, so the parameterization used here for 
that region is chosen to have a polynomial form and to 
satisfy the proper crossing symmetry while maintaining 
polynomiality for the first moments. This constrains the 
X — C dependence through the sum rule Eq[50| Having 
fixed the parameters of the Regge and diquark functions 
(Table I) , the set of measured DVCS observables can be 
determined using evolution equations to match the 
of different data sets. From these GPDs, the Compton 
Form Factors that enter the cross sections and asymme- 
tries, can be computed. Beam asymmetry data indicate 
the need for a damping of the higher C, behavior of the 
contributing CFFs. We incorporate this effect through a 
multiplicative function of X, t that lowers the high C, 
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FIG. 20: Coefficients of the beam charge asymmetry, Ac 
extracted from experiment 52, ^J'. The lower panel is the 
coefficient for the cos^ dependent term in Eq.(76c 
the upper panel is the cos independent term 
relevance of the pure DVCS contribution, | Tdvcs 
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Notice the 
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value at higher t (Eq 33 ) 



The final determination of the parameters provides an 
excellent fit to all of the available DVCS data. Newer 
DVCS cross section and asymmetry measurements at Jef- 
ferson Lab and COMPASS will provide a test of the flex- 
ibility of the model developed here. At this point we see 
that this physically motivated model provides a far reach- 
ing interpretation of the separate spin-dependent GPDs 
and thereby, a picture of the transverse structure of the 
nucleons will emerge. 

A number of questions remain that are being addressed 
in ongoing work. The connection between the Regge-like 
behavior of these GPDs and the more general form of 
variable mass diquark exchanges has opened up the possi- 
bility of having the Regge behavior emerge from diquark 
mass variation. Such a variation will better approximate 
the Fock space structure of the nucleon. A second, im- 
portant concern is the inclusion of sea quarks, whose con- 
tribution will affect the low xsj dependence, particularly 
the singlet, crossing even GPDs, whose Regge behavior is 
dominated by Pomeron exchange. Finally, the important 
extension of this parameterization scheme to the chiral 
odd GPDs is critical for the phenomenology of Deeply 
Virtual Meson Electroproduction, which was begun par- 
ticularly for the 7r° in Ref. [51] . The connection of chiral 
odd GPDs to the transversity structure of the nucleon is 



FIG. 21: Coefficients of the beam charge asymmetry, Aut, 
extracted from experiment [521 153| . The upper panel shows 
the terms E and F from Eqs.(??) and (76el, respectively; 
the middle panel shows G, and the lower panel H, both in 
Eq.(76el. The curves are predictions obtained extending our 



quantitative fit of Jefferson lab data to the Hermes set of 
observables. 



of great interest as a signal of quark and gluon orbital 
angular momentum. 
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Appendix A: Useful Integration Formulae 

A crucial point regarding the At dependence in 
Eqs.( 26|27|28|29 ) is that the integral over k_L can be done 
explicitly over the azymuthal angle first, to yield the an- 
gular dependence. This is clear in noting that fci -I- ik2 = 
Ikj^le'-^and ki-fikz = |kj^|e*'^-(l-A)/(l-C)Ai, where 
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Aj^ can be chosen to be in the x-direction with no loss 
of generahty. Also A_l -kj, = Ai|kj_| cos The 4> depen- 
dence comes only from the kj^ in the helicity flip numer- 
ators and the k'"^ = fc^ -f A^ - fc+A" fc" A+ - 2A_l • k_L 
in the denominators. When doing the integral over (p 
from to 2n the single factor e*'^gives 0, as does coscf) 
or sin(j) alone. Hence single flip amplitudes will begin at 
A terms and double flip ones at A^ . The integrals used 
are of the form 
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where b = 2A;_lAj^, and a does not depend on Aj^. 



Appendix B: Initial scale generalized parton 
distributions 



We present here a more practical version of the diquark 
contribution to the GPDs at the initial scale, = 0.0936 
GeV2. The complete parametrization needs to be multi- 



plied by the Regge term provided in Eq.(33). By defining 



L2(A:) = XMj. + {1- X)Ml- X{1- X)M^la) 
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B = 2[L^iX')-{l-X'fAl] 

all four GPDs can be written as follows, 
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For H, 
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In the forward limit the GPD's reduce to the form 



00 
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, 9O+9IK. + 92K 



where the functions 90,91,92 are evaluated at = 
and C = 0. Solving the integral for each case we get 
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Appendix C: Principal Value Integration 

The PV integrations entering the GFFs defined in 
Eq. ( 3b I were carried out with the modified Gaussian 
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method, yielding with [TU] 

J X - c J X 

C/2 C/2 

1-C } H+iX Ct)~mC Ct) ^i^"^) ^ ^^^(i/ + 2ifcos(/.) (D4) 



(CI) 



c/2 



where 



Appendix D: Kinematical coefficients in 
asymmetries and cross sections 



The kinematical coefficients in Eqs.(|77|) are 

~SK{2-y) 



1 2{2-2y + y^) 

Kdvcs = ^^^3^^4(1 -a;^,) (D2) 



^" = 3pLpUV^^? andif2^-t/Q2 
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